Resonant cyclotron acceleration of particles 
by a time periodic singular flux tube 

J. Aschi, T. Kalvoda^, P. Stovfcek^ 

3 : 1 Centre de Physique Theorique (CPT-CNRS UMR 6207) Universite du Sud, 

(N ; Toulon- Var, BP 20132, F-83957 La Garde Cedex, France (asch@cpt.univ-nirs.fr). 

Supported by Grant PHC-Barrande No. 21899QB, Ministere des affaires etrangeres 
et europeennes. 



^Department of Theoretical Computer Science, Faculty of Information Technology, 
Czech Technical University in Prague, Kolejni 2, 120 00 Praha, Czech Republic 
^ I (tomas . kalvodaOf it . cvut . cz) . 

^1 Supported by Grant No. LC06002 of the Ministry of Education of the Czech 

^ • Republic and by Grant No. 202/08/H072 of the Czech Science Foundation. 

^Department of Mathematics, Faculty of Nuclear Science, Czech Technical 
University in Prague, Trojanova 13, 120 00 Praha, Czech Republic 
^ i (stovicekOkmlinux.f jf i . cvut . cz). 

Supported by Grant No. 201/09/0811 of the Czech Science Foundation. 

^ . 

m 

^ Abstract 

We study the dynamics of a classical nonrelativistic charged particle moving 
on a punctured plane under the influence of a homogeneous magnetic field and 
driven by a periodically time-dependent singular flux tube through the hole. 
We observe an effect of resonance of the flux and cyclotron frequencies. The 
particle is accelerated to arbitrarily high energies even by a flux of small field 
strength which is not necessarily encircled by the cyclotron orbit; the cyclotron 
5^ I orbits blow up and the particle oscillates between the hole and infinity. We 

support this observation by an analytic study of an approximation for small 
amplitudes of the flux which is obtained with the aid of averaging methods. 
This way we derive asymptotic formulas that are afterwards shown to represent 
a good description of the accelerated motion even for fluxes which are not nec- 
essarily small. More precisely, we argue that the leading asymptotic terms may 
be regarded as approximate solutions of the original system in the asymptotic 
domain as the time tends to infinity. 

Keywords: electron-cyclotron resonance, singular flux tube, averaging method, 
leading asymptotic term 

AMS subject classification: 70K28, 70K65, 34E10, 34C11, 34D05 
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1 Introduction 



Consider a classical point particle of mass m and charge e moving on the punctured 
plane \ {0} in the presence of a homogeneous magnetic field of magnitude b. Sup- 
pose further that a singular magnetic flux line whose strength ^{t) is oscillating with 
frequency Q intersects the plane at the origin. The equations of motion in phase space 
P = (M^ \ {0}) X are generated by the time-dependent Hamiltonian 

Hiq,p,t) = ^ip-eAiq,t))\ with t) = f-^ + ^ V^ (1.1) 



where (g,p) G P, t G M. Here and throughout we denote x"*" = {—X2,xi) for x = 
{xi,X2) G M^. Our aim is to understand the dynamics of this system for large times. 
Of particular interest is the growth of energy as well as the drift of the guiding center. 

Our main result in the present paper is to exhibit and to prove a resonance effect 
whose origin can qualitatively be understood as follows. The Lorentz force equals 

ehq'{t)^ + eE{t) where E{t) = ---^q ■ 

27r|g|^ 

In the induction free case when $'(t) = 0, the particle moves along a circle of fixed 
center (cyclotron orbit) with the cyclotron frequency Uc = \eb\ /m and the cyclotron 
radius depending bijectively on the energy. If $'(^) is nonzero but small then the 
energy, the frequency and the guiding center of the orbit become time dependent. For 
the time derivative of the energy computed in polar coordinates (r, 6) one finds that 

fH=^^H = e q\t) . E{t) = e'm'{t) (1.2) 

and so the acceleration rate is given by 

Hir) e I r 

7 = lim = lim - / e'{t)^'{t)dt. (1.3) 

r-s-oo r 27r T-s-oo r Jq 

One may expect that 9'{t) is close to a periodic function of frequency Uc (which is the 
case when $'(t) = 0). If the frequencies Uc and VL are resonant, i.e. if there exist indices 
n and m in the support of the Fourier transforms of 9'{t) and $'(t), respectively, such 
that nujc + mVL = then one may speculate that 7 is positive. 

Our original motivation to study this problem was to understand the dynamics 
of so called quantum charge pumps [H [H] already on a classical level and to gain a 
detailed intuition on the dynamical mechanisms in the simplest case. We suggest, 
however, that particle acceleration mechanisms might be of actual interest in various 
other models, for example in interstellar physics [31 E]. We speculate that the equations 
of motion in this domain might exhibit similar ingredients as the ones studied here. 

We remark, too, that our results could be of interest in accelerator physics. While 
the betatron principle uses a linearly time dependent flux tube to accelerate particles 
on cyclotron orbits around the flux [B], the resonance effect we observe in the present 



2 



paper has the feature that acceleration can be achieved with arbitrarily small field 
strength. A second aspect is that, in contrast to the case of a linearly increasing fiux, 
cyclotron orbits which do not encircle the fiux tube are accelerated as well. In fact, 
for the linear case it was shown in |2] that outside the fiux tube one has the usual 
drift of the guiding center, without acceleration, along the field lines of the averaged 
potential. 

In the case of resonant frequencies one readily observes an accelerated motion 
numerically, see Figure [1] for an illustration. A typical resonant trajectory is a helix- 
like curve whose center goes out at the same rate as the radius grows. At the moment, 
a complete treatment of the equations of motion is out of reach. Therefore we first 
apply a resonant averaging method to derive a Hamiltonian which is formally a first 
order approximation in a small coupling constant of the fiux tube. We then analyze 
its fiow and show a certain type of asymptotic behavior at large times with the aid of 
differential topological methods. Though this approach gives only an existence result 
we use the analysis to derive the leading asymptotic terms. 

Then, in the second stage of our analysis, these asymptotic terms are used to 
formulate a simplified version of the evolution equations. We decouple the equations 
by substituting the anticipated leading asymptotic terms into the right-hand sides. 
The decoupled system admits a rigorous asymptotic analysis whose conclusions turn 
out to be fully consistent with the formulas derived by the averaging methods. This 
holds under the assumption that the singular fiux tube has a correct sign of the time 
derivative when the particle passes next to the origin. It should be emphasized that 
the singular magnetic fiux is not assumed to be very small in the analysis of the 
decoupled system. 

The paper is organized as follows. In Section [2] we summarize our main results. 
Namely, in Subsection 12. II we introduce the Hamiltonian equations of motion in action- 
angle coordinates, in Subsection 12. 21 we study the averaged dynamics resulting from the 
Poincare-von Zeipel elimination method and in Subsection 12.31 we rederive the leading 
asymptotic terms by introducing a decoupled system. Finally, in Subsection 12.41 we 
interpret the formulas derived so far in guiding center coordinates. Proofs, derivations 
and additional details are postponed to Sections [3l HI [5] and El respectively. 

2 Main results 
2.1 The Hamiltonian 

In view of the rotational symmetry of the problem we prefer to work with polar 
coordinates, q = r{cos9,sm9). Denoting by pr and pg the momenta conjugate to r 
and 9, respectively, one has 

Pr = p-q/\q\, pe = p-q^ = q^p- 
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Figure 1: The solution q{t) of the equations of motion in the plane for t G [0, 150], 
with = 27ief{Qt), f(t) = sin(t) — (1/3) cos(2t), for the values of parameters 
e = 0.35, b = l,n = l, and with the initial conditions g(0) = (1,0), g'(0) = (0, 1.617). 

Using that the vector potential A{q,t) is proportional to one finds for the Hamil- 
tonian (11. ip in polar coordinates 

The angular momentum pg is an integral of motion and thus the analysis of the 
system effectively reduces to a one-dimensional radial motion with time-dependent 
coefficients. From now on we set e = m = 1, and so the cyclotron frequency equals b. 
Put 

a{t)=pe-^m- (2.2) 

In the radial Hamiltonian one may omit the term ba{t)/2 not contributing to the 
equations of motion and thus one arrives at the expression 

^rad(r,p.,t) = | + ^ + ^. (2.3) 

For definiteness we assume that 6 > 0. 

The time-independent Hamiltonian system, with a = const, is explicitly solvable. 
In particular, one finds the corresponding action-angle coordinates (J, ip) depending 
on a as a parameter. Substituting the given function a{t) for a one gets a time- 
dependent transformation of coordinates. This is an essential step in the analysis of 
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the time-dependent Hamiltonian fl2.3p since action-angle coordinates are appropriate 
for employing averaging methods. Postponing the derivation to Section [31 here we 
give the transformation equations: 



1/2 



(2.4) 



_ ^/bl{l + \a\) cos((p) 

Pr — — ZT/2 ' 



/+^ + V/(/+|a|)sin(yp)) 



and conversely, 



Furthermore, expressing the Hamiltonian in action-angle coordinates one obtains 

(2.8) 



g,(v;,/,t) = &/-|a(t)rarctan( . , . 

Ya//+ \a{t) \ + V/sm(v9)_ 

and the corresponding Hamiltonian equations of motion take the form 

, , cos(o9)aa' 

^' = b ^ V' ^2.9) 

2y/I{I+\a\) [21 + \a\ + 2^I{I+\a\) sin((/?) j 

/ \ci\' I kl \ / N 

/' = 1 ' ' . 2.10 

2 \ 2I+\a\+2y/I{I+\a\) sin(^) J 

Using action-angle coordinates one can give a rough qualitative description of tra- 
jectories in the resonant case. From (12. 4p we see that 

= ^ + '"-) + ^ (^+ ~ sin(v5) 

where 



yf(V7TR±y7) 



are extremal points of the trajectory (see Section [3]). As formulated more precisely in 
the sequel, resonant trajectories are characterized by a linear increase of I{t) for large 
times while ip{t) ^ bt + const. Thus, as the angle if increases the radius r oscillates 
between r_ and r_|_ (though r_, r+ themselves are also time-dependent). Moreover, 
since a(t) is bounded, / — )■ oo implies r+ — )■ oo and r_(t) = 2|a(t)|/(6r_|_(t)) — )■ as 
t — 7- +00. Therefore the trajectory in the g-plane periodically returns to the origin and 
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then again escapes far away from it while its extremal distances to the origin converge 
respectively to zero and infinity. We refer again to Figured] for a typical trajectory in 
the g-plane in the case of resonant frequencies. 

Concluding this subsection let us make more precise some assumptions. For the 
sake of definiteness we shall focus on the case when po is positive and greater than 
the amplitude of $(t)/(27r) and so a{t) is an everywhere strictly positive function. 
Let us stress, however, that this restriction is not essential for the resonance effect 
as the radial Hamiltonian (12. 3p depends only on a(t)^ and thus the sign of a{t) is 
irrelevant for the motion in the radial direction. On the other hand, as discussed in 
Subsection 12.41 the sign of a{t) determines whether the orbit encircles the singular 
magnetic flux or not. 

Furthermore, the frequency > of the singular flux tube is treated as a parameter 
of the model and so we write 

$(t) = 27^e/(^]^) (2.11) 
where / is a 27r-periodic real function possibly obeying additional assumptions. Hence 

a{t)=pe-ef{nt) 

where the coupling constant e is supposed to be positive and, if desired, playing the 
role of a small parameter. 

2.2 The dynamics generated by the first order averaged Hamil- 
tonian 

In order to study occurrences of resonant behavior we apply the Poincare-von Zeipel 
elimination method which takes into account possible resonances, as explained in 
detail, for instance, in [1]. In Proposition 12 . 2 1 we provide a detailed information on the 
resonance effect for the dynamics generated by the first-order averaged Hamiltonian. 

We start from introducing some basic notation. Let T'^ = be the d- 

dimensional torus. For /(v^) G C{T'^) and A; G Z'^ we denote the /cth Fourier coefficient 
of / by the symbol 

We introduce supp as the set of indices corresponding to nonzero Fourier coeffi- 
cients of f{if). For / G CiJ'^) and L C Z'^ put 

feeL 

For example, in the formulation of Proposition 12.21 we use the averaged function 
{f{f))zu, with z/ G N. Assuming that / G C(T^) let us note that it can alternatively 
be expressed, without directly referring to the Fourier series, as 




(2.12) 
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To verify the formula observe that the RHS of fl2.12p again belongs to C(T^). To 
check that its Fourier series coincides with the LHS of fl2.12p it suffices to consider 
exponential functions f{ip) = e^^''^, with /c G Z. In that case one finds that the RHS 
of fl2.12p reproduces the function f{ip) if z/ is a divisor of k and vanishes otherwise. 

In this subsection we assume that is given by (12. lip where e > is regarded 
as a small parameter and the 27r-periodic real function f{ip) fulfills 

oo 

J2kmf{^)]k\<oo. (2.13) 

k=l 

This implies that f eC^{T^). 

Following a standard approach to time- dependent Hamiltonian systems we pass 
to the extended phase space by introducing a new phase (p2 = fit and its conjugate 
momentum I2 thus obtaining an equivalent autonomous system on a larger space. To 
have a unified notation we rename the old variables (p, I as ipi, Ii, respectively, and set 
cji = b, CU2 = ^- With this new notation, we write / = (/i, 12),^ = (v^i, ^2) (changing 
this way the meaning of the symbols / and ip in the current subsection). 

The Hamiltonian on the extended phase space is defined as 

K{ipi,ip2,Il,l2) = UJ2I2 + Hc{ipi,Il,(p2/(^2) 

Recalling (12. 8p and the conventions for pg and $ one gets 

K{y^, I, e) = uJi + UJ2I2 + eF{y^, I, e) (2.14) 

where 

F(v.,/,e) = u;,/'(^2)arctanf^=^^^4^il— ---V (2.15) 

V V ^1 +Pe-ej m) + Vh sin(v3i) J 

Our discussion focuses on the resonant case when 

A := — = — , with u, z/ G N coprime, (2.16) 

and, moreover, v fulfills 

supp ^[/] n (Zz/ \ {0}) 7^ 0. (2.17) 

Note that (I2.17P happens if and only if {f{ip))zu is not a constant function. Discussion 
of the nonresonant case is, at least on the level of the averaged dynamics, simple and 
we avoid it. Let us just remark that, in that case, it is not difficult to see that 
trajectories as well as the energy for the averaged Hamiltonian are bounded. 

The Hamiltonian (I2.14p is appropriate for application of the Poincare-von Zeipel 
elimination method whose basic scheme is briefiy recalled in Subsection 14.11 The 
idea is to eliminate from K{ip,I,e), with the aid of a canonical transformation, a 
subgroup of angle variables which are classified as nonresonant, and thus to arrive 
at a new reduced Hamiltonian IClip, J, e) depending only on the so called resonant 
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angle variables and the corresponding actions. Note, however, that a good deal of 
information about the system is contained in the canonical transformation itself. 

To achieve this goal one works with formal power series in the parameter e, and 
the construction is in fact an infinite recurrence. In particular, fC{ip, J, e) = JCoii^, J) + 
e/Ci (■?/', J) + . . . , where the leading term remains untouched, 

/Co(^, J) = M^Ij, J) = UJiJi + U2J2- (2.18) 

In practice one has to interrupt the recurrence at some order which means replacing 
the true Hamiltonian system by an approximate averaged system. In our case we shall 
be content with a truncation at the first order. 

In our model, one can apply the substitution xi = /^V'l ~ ^"02, X2 = vipi + iiip2 
where jx^v (z'L are such that ftfi + uu = 1 (/i, i> exist since fi, v are coprime). The phase 
X2 is classified as nonresonant and can be eliminated while the phase Xi is resonant 
and survives the canonical transformation. The procedure is explained in more detail 
in Subsections 14.21 and 14.31 Recalling f l2.18p combined with f l2.16p . here we just give 
the resulting formula for the averaged Hamiltonian truncated at the first order, 

/C(i)(V', J,e) = ^(i/Ji + M) + e'^i(^,^) (2.19) 

where 

/Ci(V^, J) = -^^ ^ ^[/]_„^i"^/3(Ji)l"l'^e^"(^'^^-^'^2) (2.20) 

r!,GZ\{0} 

and 



WO^,/^. (2.21) 

One observes that the averaged four- dimensional Hamiltonian system admits a 
reduction to a two-dimensional subsystem. Indeed, as it should be, }Ci{ip, J) depends 
on the angles only through the combination fl^|Jl — vil)2 = Xi- It follows that the action 
vJi + 11 J2 is an integral of motion for the Hamiltonian /C(i) (■?/', J, e). The reduced 
system then depends on the coordinates Xi, Ji- 

Moreover, by inspection of the series fl2.20p one finds that the Hamiltonian for 
the reduced subsystem can be expressed in terms of a single complex variable z = 
/3( Ji)'^e*-^^ and takes the form Z{xi,Ji) = Re[/;,(2;)], with h being a holomorphic 
function. We refer to Subsection 14.31 for the proof of the following abstract theorem. 

Theorem 2.1. Lethe C^(BI) and suppose h{z) is a nonconstant holomorphic func- 
tion on the open unit disk Bi. Let g : [0, +oo[ — )■ [0, 1[ &e a smooth function such that 
q'{x) > for X > 0, g{0) = and lima;._>.+oo q{x) = 1. Let Z{xiy Ji) be the Hamilton 
function on Mx ]0, +oo[ defined by 

Z(xi,Ji) = Re[/i(f?(Ji)e*>^0]- 

Then, for almost all initial conditions (xi(0), Ji(0)), the corresponding Hamiltonian 
trajectory fulfills 

lim xi{t) = Xi(oo) e lim Ji{t) = +00, (2.22) 
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and 



lim J[{t) = Im[e*^^(~) /i'(e^^^(°°))] > 0. 



(2.23) 



Theorem 12.11 yields the desired information about the asymptotic behavior of the 
averaged Hamiltonian system in action-angle coordinates {J,ip)- One has to go back 
to the original action-angle coordinates (/, (p), however, since in these coordinates the 
dynamics of the studied system can be directly interpreted. This means to invert the 
canonical transformation (/, v^) {J,ip) which resulted from the Poincare-von Zeipel 
elimination method. Let us remark that the generating function of this canonical 
transformation is truncated at the first order as well. Doing so one derives the following 
result whose proof can again be found in Subsection 14.31 

Proposition 2.2. Suppose that conditions ^2.16\) . 1^2. 17\ ) characterizing the resonant 
case are satisfied. Let {ipi{t), hit)) be a trajectory for the first-order averaged Hamil- 
tonian i\2.iy\) but expressed in the original angle-action coordinates, as introduced in 
Subsection \2.1\ Then, for almost all initial conditions (v9i(0), Ji(0)), 



2.3 A decoupled system and the leading asymptotic terms 

In the preceding subsection we studied an approximate dynamics derived with the 
aid of averaging methods. It turned out that a typical trajectory reaches for large 
times an asymptotic domain characterized by formulas 02.241) . Being guided by this 
experience as well as by numerical experiments that we have carried out we suggest 
that formulas fl2.24p in fact give the correct leading asymptotic terms for the complete 
system. 

To support this suggestion we now derive the leading asymptotic terms and a 
bound on the order of error terms for the original Hamiltonian equations (12. 9p . (I2.10p 
while assuming that the system already follows the anticipated asymptotic regime. At 
the same time, we relax the assumption on smallness of the parameter e. 

Applying the substitution 




C > 0, 



(2.24) 



and 



c = --^ /:(- (0(00) + f ) a) , with /.(y,) = (/(v,)) 



F{t) = 21{t) + |a(t)|, 0(t) = ifit) - bt, 



to equations (12. 9p . (I2.10p we obtain the system of differential equations 



F'{t) 



a{t) a' it) 



(2.25) 



0'(t) 



F{t) + ^F(t)2 -a(t)2sin(6t + 0(t)) ' 

a{t) a'{t) cos(6t + 0(t)) 



(2.26) 



v/F(t)2 - a(t)2 (F{t) + v/F(t)2-a(t)2sin(6t + 0(t)) 
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where a{t) = pe — €f{Qt) (see (12. 2p and (I2.1ip ). The real function f{t) is supposed 
to be continuously differentiable and 27r-periodic. The only assumptions imposed on 
e > in this subsection are that e is not too big so that the function a{t) has no zeroes 
and is everywhere of the same sign as pe. Recall that for definiteness pe is supposed 
to be positive. Clearly, the functions a{t) and a'{t) are bounded on M. 

Equations (I2.25P and (I2.26P are nonlinear and coupled together. To decouple 
them we replace 0(t) on the RHS of dMSD and F{t) on the RHS of (^M) by the 
respective leading asymptotic terms (I2.24p . as learned from the averaging method. 
This is done under the assumption that the solution has already reached the domain 
F{t) > Fq ^ pg where F[t) is sufficiently large and starts to grow. Thus, to formulate 
a problem with decoupled equations, we replace (j){t) in (I2.25p by the expected limit 
value (p = 0(00) G M, i.e. the simplified equation reads 

^// N ciit) a' it) 

F'it) = ^ ^ ^ ^ . 2.27 

F{t) + ^/F{ty - a(t)2 sm{bt + 0) 

As stated in Proposition 12.31 below, a solution F{t) of (12.271) actually grows linearly 
for large times. Conversely, equation (I2.26P is analyzed in Proposition 12.41 under the 
assumption that F{t) grows linearly. In that case (j){t) is actually shown to approach a 
constant value as t tends to infinity. For derivation of Proposition 12.41 the periodicity 
of the functions a{t) is not important. It suffices to assume that it takes values from 
a bounded interval separated from zero. For the proofs see Section [51 

Proposition 2.3. Suppose f G C^(M) is 2Tt -periodic and 



Then any solution of l^2.21\j such that -F(O) is sufficiently large fulfills 



F{t) = en 



t + 0{\og{tf) as t^ +00. 



Proposition 2.4. Suppose a{t) E C^(R) fulfills 

Ai <a{t) < A2, |a(t)| </l3, 

for some positive constants Ai, A2, A3. Furthermore, suppose F{t) G C(M) has the 
asymptotic behavior 

F{t) = at + 0{log{tf) as t-^ +00, (2.28) 

with a positive constant a. If (pit) obeys the differential equation 26]) on a neigh- 
borhood 0/+00, then there exists a finite limit limt^^oo 4>{t) = 4>{oo) and 

(/)(i) = 0(00) + O ( ) as t^ +00. (2.29) 
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Remarks, (i) Let us point out an essential difference between equations fl2.25l) 
and fl2.26p . Note that for all F, a, s G M such that F > \a\ > one has 



acos s 



F + y/F^-a^ sin(s) 
and, consequently, it follows from f l2.26p that 



< 1, (2.30) 



Thus the RHS of f l2.26p is inversely proportional to F(t). The point is that if the 
denominator in f l2.30p becomes very small for those s for which sin(s) = —1 then one 
gets a compensation by the vanishing numerator. Nothing similar can be claimed, 
however, for equation fl2.25p . 

(ii) The replacement of the phase (f){t) by a constant in fl2.25p was quite crucial for 
derivation of the result stated in Proposition 12. 3[ In fact, suppose that F(t) is suffi- 
ciently large. In the case of the original equation f l2.25p . too, essential contributions to 
the increase of F{t) are achieved at the moments of time for which sm{bt + (j){t)) = —1. 
If (j){t) equals a constant then these moments of time are well defined and the growth 
of F{t) can be estimated. On the contrary, without a sufficiently precise information 
about one loses any control on the growth of F{t). 

Propositions 12.31 and 12.41 can also be interpreted in the following way. Let us pass 
from the differential equations (12.25 p . f l2.26p to the integral equations 

/"* a(s) a'(s) 
Ft -FO - / ^ ' ^ ds = 0, 

- 0(oo) - r a{s)a'{s)c os{bs + <P{s)) ^ ^ ^ 

Jt ^/F{sy - a(s)2 (f{s) + ^F{sy - a{sY sin(6s + 

Suppose 0(oo) satisfies a = -e fi/'(-(fi/6)(0(cx)) + n/2)) > 0. If F(0) is sufficiently 
large then the functions 

F{t) = at + F(0), (j){t) = 0(oo), t > 0, 

can be regarded as an approximate solution of this system of integral equations with 
errors of order 0(log(t)2) for the first equation and of order 0(log(t)/t) for the second 
one. 

One has to admit, however, that this argument still does not represent a com- 
plete mathematical proof of the asymptotic behavior of the action-angle variables 
I{t) = (F(t) — |a(t)|)/2, (p(t) = bt + (f){t). So far we have analyzed either the dynamics 
generated by an approximate averaged Hamiltonian in Subsection 12.21 or a simplified 
decoupled system in the current subsection. Moreover, in the latter case it should be 
emphasized that the simplified equations were derived under the essential assumption 
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that the dynamical system had already reached the regime characterized by an ac- 
celeration with an unlimited energy growth (this is reflected by the assumption that 
F(0) is sufficiently large). Nevertheless, on the basis of this analysis as well as on the 
basis of numerical experiments we formulate the following conjecture. 

Conjecture 2.5. If ^l/b G N then the regime of acceleration for the original (true) 
dynamical system is described by the asymptotic behavior 

i{t) = ct + o (iog(t)2) , ip{t) = bt + (j){oo) + o (^^^Y^^ ^ ^ ^^-^^^ 

where 0(oo) is a real constant and 

C = -\ eilfi-O > 0, wzth e = J (0(00) + I) . (2.33) 

2.4 Guiding center coordinates 

Being given the asymptotic relations fl2.32p . f l2.33p it is desirable to describe the ac- 
celerated motion in terms of the original Cartesian coordinates q. The description 
becomes more transparent if the motion is decomposed into a motion of the guiding 
center X and a relative motion of the particle with respect to this center which is 
characterized by a gyroradius vector R and a gyrophase |[8j. Since the presented 
results have a direct physical interpretation, in this subsection we make an exception 
and give the formulas using now all physical constants (including m and e). Additional 
details and derivations are postponed to Section O 

Let v{q,p,t) = {p — eA{q,t))/m be the velocity. We write q = X + R where, by 
definition, 

X{q,p,t) = q + — v^{q,p,t), R{q,p,t) = q - X{q,p,t) = — -v^{q,p,t) 

Uc UJc 

are the guiding center field and the gyroradius field, respectively. In what follows, we 
use the polar decompositions 

g = r (cos 6^, sin 6^), X = |X|(cosx, sinx), i? = ji^Kcos-i?, sint?). 

Concerning the geometrical arrangement, one has the relation 

■^ = ip + ^-l (mod 27r) (2.34) 

(with (f being introduced in Subsection 12. 1^ for the derivation see Section E]). 

The quantities X, R were introduced (under different names) and studied in [2] 
where one can also find several formulas given below, notably those given in f l2.36p . 
Observe that |pp = p^ + Po/r'^ and, according to (12.31) and (12. 6p . 

H = H,^A + — = ujJl + -i\a\+a)] , with a{t) = pg <^{t). (2.35) 



2 V. . . , > ^ ' ^ " 27r 
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Using these relations one derives that 

!^^|xP = i/-c..p„ + ^^*, = (2.36) 

Observe from (12.361) that if a{t) is an everywhere positive function then \R{t)\ > \X{t)\ 
and so the center of coordinates always stays in the domain encircled by the spiral-like 
trajectory. On the contrary, if a{t) is an everywhere negative function then the center 
of coordinates is never encircled by the trajectory. 

For a Hamiltonian trajectory {q{t),p(t)) put H(t) = H{q(t),p(t),t). Suppose again 
that pe > and hence a(t) > 0. Applying Conjecture 12.51 and recaUing (12.111) one has 

/(t) = -^$'(^--i(^0(oo) + ^)^t + O(log(t)2) ast^+oo. (2.37) 

Using (I2.35P and recalling definition (II. 3p of the acceleration rate 7 one finds the 
positive value 

From (12. 36 p . (I2.35P and (I2.37P one also deduces the asymptotic behavior of the guiding 
center and the gyroradius, 



|X(t)| = J^+0(M2!Y|flM| = j^+0(M2!'l (2.39) 



Now it is obvious why the resonance has also consequences for the drift. The growing 
energy H comes along with a growing distance of the guiding center from the origin. 
Thus in g-space the particle oscillates around a drifting center which goes out at the 
same rate as the radius grows. 

From Conjecture 12. 5[ with a bit of additional analysis (see Section |6]), one can also 
derive consequences for the asymptotic behavior of the angle variables. One has, as 
t — > +00, 

x{t) = D\og{u,t)+x{oo) + o{l), m =u,t + D log{uj,t)+ ^{00) + o{l), (2.40) 

where D, x(oo) and ^{00) are real constants. Relations (12.391) and (I2.40p give a 
complete information about the asymptotic behavior of the trajectory q(t) = X{t) + 
R{t). 

Finally, in connection with formula (I2.38P let us point out that the studied system 
behaves for large times almost as a "kicked" system. Note that, for a fixed a > 0, 

\n\/2 



X — y/x'^ — a'^ cos(t) 



> e*"* 27r527r(i) as a; +00, 

^-^ \x + a J 
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where S^it) = Yln=-oo^(^ ~ stands for the A-periodic prolongation of the Dirac 6 
function. Since H'{t) = -e^'(t)<l>'(t)/(27r) (see ([12])), ^(t) ^ + 0(oo) and, in view 
of (J22D, 



one has 




Thus the system gains the main contributions to the energy growth in narrow intervals 
around the instances of time for which uj. = — </)(oo) — 7r/2 (mod 27r). According to 
f l2.34p . these are exactly those instances of time when ■(9(t) — x{t) ~ (mod 2tt), i.e. 
when the particle passes next to the singular flux line. 



3 Derivation of the Hamiltonian in action-angle co- 
ordinates 

In this section we derive the transformation equations from coordinates r, to action- 
angle coordinates I, (p. 

Let us first note that the equations of motion for the radial Hamiltonian (12. 3p have 
a well-defined unique solution for any initial condition and for all times. They are 
equivalent to the nonlinear second-order differential equation 



62 



a{ty 



(3.1) 



Proposition 3.1. Suppose a{t) is a real continuously differentiable function defined 
on M having no zeros. Then for any initial condition r{to) = tq, r'(to) = ^i, with 
(to,ro,ri) G Mx ]0, -|-oo[ xM, there exists a unique solution of the differential equation 
lis. 1\} defined on the whole real line M and satisfying this initial condition. 

Put pr = r' and 



Proof. Suppose r(t) is a solution of the Hamiltonian equations (13.1 
H{t)=H,^{r{t),pr{t),t). Then 



\a{t)a'{t)\ 
r(t)2 



< 



2a'{t) 



a{t) 



Hit). 



From here one readily concludes that if r{t) is defined on a bounded interval M C M, 
then there exist constants Ri, R2, < Ri < R2 < +00, such that Ri < r(t) < R2 for 
all t G M. From the general theory of ordinary differential equations it follows that 
any solution r{t) of (13. ip can be continued to the whole real line. □ 

As a first step of the derivation of the transformation we introduce the action-angle 
coordinates for a frozen time. Assume for a moment that a{t) = a is a constant and 
denote 



,2^2 



V{r) 



2r^ 



+ 



b''r 



14 



Suppose a fixed energy level E is greater than the minimal value Knin = b\a\/2. Then 
the motion is constrained to a bounded interval [r_,r_|_], and one has 

E - Vir) = (r 2 - r2)(r2 - r^). 



8r2 



The action equals 



I{E) = - r ^2{E-V{p))dp 

71" Jr^ 



r 1 — r_ 



Hence 



Using the generating function 



1/2 



yf(v/7TR±v/7). 



5(r 



one derives the canonical transformation of variables between (r,pr) and the angle- 
action variables {^p,I)- One has 



as 
a7 



6 
4 

TT 



^ 2 _|_ ^ 2 _ 



p2 _ ^2 



dp 



arctan 



■Z (2^/^ -r2)(r2 -r2) 

For the angle variable = dS/dl — 7t/2 one obtains 



sin((y9) 



VHl+\a\) V 4 



Furthermore, 



^ / 2 



Relations (El]), (El]), ([M]), ([22]) readily follow. 

Let us now switch to the time-dependent case. The Hamiltonian transforms ac- 
cording to the rule 



dS{u,I,t) 



dt 



u=r{ip,I ,t) 



One computes 
dS{u,I,t) 



dt 



b \a\' r(^'^'*) 1 / /p2 _ ^2 1,2 _ p2 



u=riip,I,t) 



p rf — p'^ d\a\ y p^ — r"^ d\a\ 



dp. 



Simplifying the expression and dropping those terms which are independent of if and 
/ one finally arrives at the expression (12.81) . 
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4 Details on the von Zeipel elimination method 
and the proofs 

4.1 A summary of basic formulas 

This is a short summary of basic steps in the Poincare-von Zeipel ehmination proce- 
dure. Recall the notation from the beginning of Subsection 12.21 

Consider a completely integrable Hamiltonian in action-angle coordinates, Kq{I) = 
u ■ L where / runs over a domain in M'^, If e T'^ and u G is a constant vector 
of frequencies. One is interested in a perturbed system with a small Hamiltonian 
perturbation so that the total Hamiltonian reads 

K{ip, I, e) = Ko{I) + e K,{^, I, e) = Ko{I) + e K^{ip, I) + K,{^, I) + . . . 

where e is a small parameter. The function K^[ip, I, e) is assumed to be analytic in all 
variables 

Let K be the lattice of indices in Z"^ corresponding to resonant frequencies and W 
be its complement, i.e. 

K= {n eZ"^; n-u = 0}, W = Z'^ \ K. 

One applies a formal canonical transformation of variables, {I,<f) H- {J tip)-, so that 
the Fourier series in the angle variables ip of the resulting Hamiltonian }C{ip, J, e) has 
nonzero coefficients only for indices from the lattice K. The canonical transformation 
is generated by a function S{(f, J, e) regarded as a formal power series with coefficient 
functions Sj{ip, J) and the absolute term So{(p, J) = ip- J. Similarly, the new Hamilto- 
nian }C{tp, J, e) is sought in the form of a formal power series with coefficient functions 
}Cj{ip, J). One arrives at the system of equations /Co(J, (p) = Ko{J) = u ■ J and 

]CM J) = ^- ^^^^ + Jl J > h 

where Pi{ip, J) = Ki{(p, J) and the terms Pj for j >2 are determined recursively. The 
formal von Zeipel Hamiltonian is defined by the equalities JCj{ip, J) = {Pji'ip, </))k for 
j > 1. Coefficients Sj{ip, J) are then solutions of the first order differential equations 

u.-^^^ = -{PMJ))^.,J>l. 

In practice one truncates /C ('?/', J, e) at some order m > 1 of the parameter e. Let 
us define the mth order averaged Hamiltonian 

/C(^)(^, J, e) = /Co(J) + e/Ci(^, J) + . . . + e'"/C„(V', J). 

Similarly, let S'(m)((y9, J, e) be the truncated generating function. If {ip(t),J{t)) is a 
solution of the Hamiltonian equations for /C(m,) ('?/', J, e), and if {(fit), I(t)) is the same 
solution after the inverted canonical transformation generated by S'(m)(v^, ^, e), then 
{(p{t),I{t)) is expected to approximate well the solution of the original system (gov- 
erned by the Hamiltonian K{ip,I,e)) for times of order l/e™ (see [1] for a detailed 
discussion) . 
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4.2 Derivation of the first-order averaged Hamiltonian 

If the ratio UJ2/UJ1 is irrational then the lattice IK is trivial, IK = {0}, and the von 
Zeipel method amounts to the ordinary averaging method in the angle variables if. 
Here we focus on the complementary case fl2.16p when bj^ji^x = /i/'^, with /i and u 
being coprime positive integers. 

Recalling din]), we have /^(v?, /, e) = U1I1 + UJ2I2 + eKi{ip, I) +e'^K{ip, I ,e) 

where K{ip, I, e) is an analytic function in e, 

K,{^,I) = uj2nv2)F,{viJi), = arctan f-==^^^^^p^-—-y 

Vv/i +pe + \/h sm(v9i)/ 



One finds that 



^[Fi(y.i,/i)]fc = ^(^-^-^J forA;^0, /i)]o = 0. 

Obviously, the Fourier image of Ki{ip,I) takes nonzero values only for indices {k,i) 
such that keZ\ {0}, i G supp ^[f] \ {0}, and 

Now we can describe the von Zeipel canonical transformation of the first order. 
The resonant lattice is IK = Z (/i, — z^), and, according to the general scheme, one has 
]Ci{ip,J) = {Ki^ip, J))^. With the above computations, this immediately leads to 
(I22DD, fl^:^ and, consequently, to (KT^ . 

An important piece of information is also the Poincare-von Zeipel canonical trans- 
formation generated by S(^i){ip, J, e). Si{ip, J) is a solution to the differential equation 
CO ■ dSi/dip = —Ki + /Ci. Seeking Si{ip, J) in the form 

Si{^, ^) = 2 Re 1^ '^[f'VGkivi. Ji) e''^^ j (4.1) 
one finally arrives at the countable system of equations 

^ ^'^ ^ neZ\{0}, n^-kX 

For the solution we choose 

neZ\{0}, n=/=-k\ ^ ' 

Of course, \i k\ ^ 'L then the restriction n ^ —kX is void. On the other hand, if 
k\ G Z, and this happens if and only if A; G Zz/, then the solution Gk{^Pi, Ji) is not 
unique. 
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4.3 Proofs of Theorem 12.11 and Proposition 12.21 

The purpose of this subsection is to give some details concerning the dynamics gen- 
erated by the Hamiltonian ]C(^i-){ip, J, e), as defined in fl2.19p . and the proofs of Theo- 
rem 12.11 and Proposition 12.21 

We start from the reduction to a two-dimensional Hamiltonian subsystem. Since 
fi and z/ are coprime there exist /i, z/ G Z such that p.^ + uu = 1. Put 



R 



fl —V 



and consider the canonical transformation x = R-V') J = ^J^L generated by the 
function L ■ R?/;. In particular, 

Since }Ci{4', J) in fl2.20p depends only on the angle xi; and not on X2, the action L2 
is an integral of motion. Let us define 

Then 

9/Ci(^, J) dJ, dZ{xi,Ji) 



m = 

J[{t) = -e 



dJi dLi dJi 
dJC^i^P, J) 1 dZixu Jl) dxi dZixu Jl] 



d^Ji /i dxi d^pi dxi 



Thus the time evolution in coordinates xi, Ji is governed by the Hamiltonian Z{xi, Ji)- 
From the form of of the series (12.201) one can see that the Hamiltonian Z{xi, Ji) 
can be expressed, in a compact way, in terms of a holomorphic function. Set 

00 

h{z) = -efiuiJ2'^ifU.t''^z" (4.3) 

n=l 

and 



q{x) = I3{xr = — , X > 0. 

\x + pej 

Then, by assumption f l2.13p . h{z) is holomorphic on the open unit disk i?i C C 
and h G C^{Bi). Moreover, one can write Z{xiiJi) = ^^[h{Q{Ji)e'^^^)] (note that 
'^[f]nu = ■^[f]-nu since / is real). The Hamiltonian equations of motion read 

X'^it) = 4^ Re[zh'{z)], J[{t) = lm[zh'{z)], with z = Q{Ji)e'^\ (4.4) 
Q\Ji) 

Concerning the asymptotic behavior of Hamiltonian trajectories (xi(t), Jiif)), as t — )• 
-|-oo, one can formulate a proposition under somewhat more general circumstances, as 
it is done in Theorem 12.11 
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Proof of Theorem \2.1\ Set R{z) = Iie[h{z)], z E Bi. Then 



dR, = {Re[h'{z)],-lm[h'{z)]). 

Hence dRz = iff h'{z) = 0, and the set of critical points of R in Bi has no accumu- 
lation points in Bi and is at most countable. By Sard's theorem, almost all y G M 
are regular values of R\dBi. If y is a regular value both of R and R\dBi then the 
level set R~^{y) is a compact one-dimensional submanifold with boundary in Bi, 
dR~^{y) = R~^{y) fl dBi and R~^{y) is not tangent to dBi at any point. Moreover, 
R~^{y) n i?i is a smooth submanifold of i?i [1]. By the classification of compact 
connected one-dimensional manifolds [1], every component of R~^{y) is diffeomorphic 
either to a circle or to a closed interval. But the first possibility is excluded because 
R{z) is a harmonic function. In fact, if U G Bi, U is an open set, dU ~ S"^ is a 
smooth submanifold of Bi and R{z) is constant on dU then R{z) is constant on U 
and so is h{z). Consequently, h{z) is constant on Si, a contradiction with our as- 
sumptions. Thus every component F of R~^{y) is diffeomorphic to a closed interval, 
dr = {a,b} = r n dBi, and F is tangent to dBi neither at a nor at b. Let z G Bi 
be such that dRz 7^ 0. By the local submersion theorem [1], i? is locally equivalent at 
z to the canonical submersion — )■ M. Hence z possesses an open neighborhood U 
such that R{U) is an open interval. We know that almost every y G R{U) is a regular 
value both of R and R\dBi. By the Fubini theorem, for almost every w E U, R{w) is 
a regular value both of R and R\dBi. The same claim is true for almost all w E Bi 
because the set of critical points of R in Bi is at most countable. It follows that for 
almost all {xi, Ji) ^ Mx]0,+oo[, R{g{Ji)e^^^) 7^ R{0) is a regular value both of R 
and R\dBi. 

Suppose now that an initial condition (xi(0), -'^i(O)) has been chosen so that y = 
i?(^(Ji(0))e*'^i'^°)) 7^ -R(O) is a regular value both of R and R\dBi. Let F be the compo- 
nent of R~^{y) containing the point g{Ji{0))e^^^^^\ Since the Hamiltonian Z(xi, Ji) is 
an integral of motion the Hamiltonian trajectory z{t) = Ji(t))e*^^''*^ is constrained 
to the submanifold F C -Bi. We have to show that z{t) reaches the boundary dBi as 
t — )■ +CO. The tangent vector to the trajectory at the point z{t) equals 



-^ = iQ{Ji{t))Q\Ut))h'{z{t)). 

Since ^ F, q\Ji) > for all Ji > and h'{z) has no zeroes on F (because y is a 
regular value) it follows that z{t) leaves any compact subset of Bi in a finite time. It 
remains to show that z{t) does not reach dBi in a finite time. But by equations of 
motion fl4.4p . | < max^g^^^ ^ind so Ji{t) cannot grow faster than linearly. 

This reasoning shows ([222]). From (|0]) and (12^^ it follows (OSD; one has only 
to justify the sign of the limit. Obviously, the limit must be nonnegative. Denote 
dR = R\dBi. Then dR can be regarded as a function of the angle variable, dR{x) = 
Re[/i(e*^)], and one has 

(9i?)'(xi(oo)) = -Im[e^>^i(°^)/i'(e^»(°°))] ^0 

because y = dR{xi{^)) is a regular value of dR. 



19 



Finally, to derive Proposition 12.21 one has to apply the inverted canonical transfor- 
mation, from (ip, J) to {ip, I), generated by 5'(i)(v9, J, e) = (p ■ J + eSi{(p, J). Hence 

^ = ^ + e^j—, I = J + ^^^- (4.5) 

For the proof we need a couple of auxiliary estimates. 
Lemma 4.1. For all /S, < l3 < 1, and all aeR\Z, 

E f 1 ^ A^Tf-^ + 2 + 6| log(l - (4.6) 

jez ^ \n-a\ dist(a,Z) 



and for all a G Z, 



f .<l + 3|log(l-/3)|. (4.7) 

\n — a\ 



nGZ, n^a 



Proof. Notice that inequality (14. 6 p is invariant if a is replaced either by —a or 
by A; + a. A; G Z. Thus we can restrict ourselves to the interval < a < 1/2. Then 
\a\ = dist(a, Z) and |r;, — a| > \n\/2. This observation reduces (14. 6 p to (14. 7p with a = 0. 
Similarly, inequality (14. 7p is invariant if a is replaced by + a. A; G Z. It follows that, 
in both cases, it suffices to show (14. 7p for a = and with the range of j restricted to 
nonnegative integers. 

Suppose that j > 0. Splitting the range of summation in n into the subranges 
n < — j — 1, — j < n < — 1 and 1 < n, one gets 

a\n+j\ ^ Rj—'m 

y: ^<2iiog(i-^)i+^^. 

^-^ \n\ ^-^ m 

aj—m am 

E E -<iiog(i-/3)i 

^<m<j/2 l<m<j/2 

m ^-^ m 

j/2<m<j j/2<m<j 

This shows the lemma. 

Proof of Proposition \2.2[ If (12.171) is true then h{z) defined in (14.31) obeys the 
assumptions of Theorem 12.11 and so for almost all initial conditions (xi(0), Ji(0)), 
equalities (I2.22p and (I2.23p hold. In particular. Theorem 12. II implies that 

1 -/3(Ji(t)) = O(r^) as t ^ +00. (4.8) 

Putting f^{ip) = {f{f))zu one also has 

hm = Imre^^^(-) /.'(e^^^(-))l = fJ-^^ -—)>0. (4.9) 



Furthermore, 



and 
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Further, using fl4.2p one can estimate (recalling that A = ^i/v) 



dJi 

Hence, using (12.211) . 



raeZ\{0}, n=/=-k\ 



n + kX 



< 



l-/3(Ji) 



dJi 



< 



1^ l-/3(Ji) 
2pe P{Ji) 



(4.10) 



As a next step let us estimate \dGk/ dipi\. Clearly, 



- 2 



aeZ\{0}, nj^-k\ 



\n + kX\ 



Writing k\ = —j — a, with j G Z and a = s/u, s = 0, 1, . . . , z/ — 1, one can apply 
Lemma 14.11 to show that 



dip I 



<c' + c"|log(l-/3(Ji)) 



(4.11) 



where the constants c', c" do not depend on k and (pi, Ji. 

From definitions (14. ip . (14. 2 p and from assumption (I2.13P one can readily see that 
S{ip, J) is in and C°° in ipi, Ji (and does not depend on J2). Moreover, from 
dim and (USD it follows that 



dS,{ip{t),J{t)) 
dJi 

Similarly, estimate (14. lip implies 

dSi{^{t),J{t)) 
dip I 



0{t ^) as t +00. 



0(log(t)) as t +00. 



(4.12) 



(4.13) 



Now one can deduce the asymptotic behavior of ipi{t) and Ii{t). Observe that 
11)2 = UJ2 and Si{ip, J) does not depend on J2, hence '?/'2(t) = p>2{t) = ^2t- Furthermore, 
i>i = {Xi + i^i'2)/fJ' and so 



lim (V^i(t) - uJit) = - Xi(oo). 



(4.14) 



Put 0(cxd) = xi(oo)//i. To conclude the proof it suffices to recall the transformation 
rules (14. 5 p and to take into account (14.91) jointly with (14.131) and (I4.14p jointly with 
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5 Proofs of Proposition 



2.3 



and Proposition 
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Analyzing equation f l2.27p we prefer to work with a rescaled time (or one can choose 
the units so that 6 = 1) and, simphfying the notation, we consider a differential 
equation of the form 

9'it) = (5.1) 



where g{t), a{t) are continuously differentiable real functions, a{t) is strictly positive 
and is a real constant. In the resonant case the functions g{t), a{t) are supposed to 
be 27r-periodic which means for the original data that G N (and 6 = 1). 

In the first step we estimate the growth of a solution on an interval of length 7r/2. 
Let 11/11 = max|/(t)| denote the norm in C([0,7r/2]), and put 

A = min a(t) > 0. 

0<t<7r/2 

Consider for a moment the differential equation 

h'{t) = , ^^^^ (5.2) 

h{t) - ^Jhitf-aitf cos(t) 

on the interval [0, 7r/2] with an initial condition h{{]) = Hq > \\a\\. The goal is to show 
that for large values of h^, an essential contribution to the growth of a solution h{t) 
on this interval comes from a narrow neighborhood of the point t = 0. 

Remark 5.1. If g{t) is nonnegative on the interval [0,7r/2] then a solution h{t) to 
li5.S\) surely exists and is unique. In the general case, the existence and uniqueness 
is guaranteed provided the initial condition ho is sufficiently large. From liS.^) one 
derives that \h'{t)\ < 2\\g\\h(t)/A'^ and so 

exp{-2\\Q\\t/A^) ho < h{t) < exp{2\\Q\\t/A'^) ho (5.3) 

as long as h{t) makes sense. Consequently, a sufficient condition for existence of a 
solution is ho > exp(7r||^||/y4^)||a|| . 

Lemma 5.2. Let g,a E C"'^([0, 7r/2]) be real functions, g{0) ^ and a{t) > on 
[0,7r/2]. Consider the set of solutions h{t) to the differential equation h5.2^) on the 
interval [0,7r/2] with a variable initial condition h{0) = ho for ho sufficiently large. 
Then 

h (I) =ho+ + O ( log(/io)) as ho ^ +oo. 

Proof. Let us fix < < 7r/2, so that \g{t)\ > on the interval [0,?7[, i.e. 
g{t) does not change its sign on that interval. Thus any solution h{t) to (15. 2 p is 
strictly monotone on [0,?]]. For r] < t < tt/2 one can estimate |/i'(t)| < C/h{t) where 
C = \\g\\/{l — cos{r])). In view of (15. 3p it follows that 

h{7r/2) - h{r]) = O (V^) as ho +oo. (5.4) 
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Set 

hi = m.m{h{0), h{ri)}, h2 = max{/i(0), h{ri)}, A = h{ri) — h^. 
Then |A| = /i2 — hi. Set, for x > a > 0, 

1 



X — \/ x^ — cos(t) 

One has h!{t) = g{t) a(t), t). If x > 2u/V3 > then Vx^ - > x/2 and 



I v^a;2 — — a; cos(t) I 



\/ x'^ — {x — \/. 



x^ — cos(t))" ^ 



9m 



'^{x,u,t) 



u cos(t) 



— u- 



{x — \/ x^ — v? cos(t)) 



2 < -^(x,M,t), (5.5) 
3 

< -^(x,M,t). (5.6) 



Observe also that, for x > u > 0, 



V2 2 

\E'(x, M, t) dt = — arctan 
w V 



/x + a/ 



X^ — n?' \ TT 

< - 



(5.7) 



Assuming that /iq is sufficiently large and using (15. 5p . (15. 7p one can estimate 



A - /%(t)^(/io,a(t),t)dt 





< / |f?(t)||^(/i(t),a(t),t)-^(/io,a(t),t)|dt 
'o 

< 3^ j^^^' (^j\^^^A,t)dx^dt 



In view of (15. 3p it follows that 

A= (1 + 0(V')) r gmiho,a{t),t)dt. 
Jo 

Furthermore, with the aid of (15. 6 p one finds that 

\g{t)^{ho,a{t),t) - g{0)^{ho,a{0),t)\ < C'^{ho,A,t)t 

where 



C'= 1 + 



9||g|P 



2\ 1/2 



Note that 



TT 



A, t) t dt < - 



tt/2 



sin(t) 



2 7o /^o - V/io - cos(t) 



dt = o(v'iog(M) 
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and 



^%(0)^(/io,a(0),t)dt = g{0) J ^ ^{ho,a{0),t)dt + 0{ho') 

vrp(O) 



a(0) 



Altogether this means that 



h{v) -ho = A = ^ + 0(V' log(/^o)) . 

Recalhng fl5.4p . the lemma follows. 

Consider the mapping T-t : h{0) i— )■ /i(7r/2), where runs over solutions to the 
differential equation (15. 2p . From the general theory of ordinary differential equations it 
is known that "H is a mapping well defined on a neighborhood of +00. Lemma 15^ 
claims that 'H(x) = x + 7r^(0)/a(0) + 0(x~^ log(x)). On the basis of similar argu- 
ments, the inverse mapping "H^^ : hiji/2) 1— )■ /i(0) is also well defined and on a 
neighborhood of +00. From the asymptotic behavior of 'H(x) one readily derives that 
^~^(?/) = y ~ vr^(0)/a(0) + 0{y~^\og{y)). These considerations make it possible to 
reverse the roles of the boundary points and 7r/2. Moreover, splitting the interval 
[0, 27r] into four subintervals of length 7i/2 one arrives at the following lemma. 

Lemma 5.3. Let g,a E C^([0,27r]) be real functions, g{7r) 7^ and a{t) > on [0,27r]. 
Consider the set of solutions h{t) to the differential equation 

h'it) '^'^ 



h{t) + ^h{t)^ - a(t)2 cos(t) 

on the interval [0, 2tt] with a variable initial condition h{0) = Hq for Hq sufficiently 
large. Then 

h{27T) = ho + + o ( log(/io)) as ho +00. 

a[TT) 

In the next step, applying repeatedly Lemma 15.31 one can show that solutions of 
the differential equation (15.1 p in the resonant case f2 G N (with 6 = 1) grow with time 
linearly provided the initial condition is sufficiently large and the phase belongs to 
a certain interval. 

Proposition 5.4. Suppose g{t), a{t) are continuously differentiable 27^ -periodic real 
functions, G M, a{t) is everywhere positive and 
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Let g{t) be a solution of the differential equation / 1 5. 1\) on the interval t > with the 
initial condition g{0) = go>l- If go is sufficiently large then 

Proof of Proposition \2.3[ This is just an immediate corollary of Proposition 15.41 
It suffices to go back to the original notation and equation fl2.27p by applying the 
substitution F{t) = g{ht), a{t) = d{bt), g(t) = a{t)a'{t). Equation fl2.27p transforms 
into fl5.ll) (with a(t) being replaced by a(t)), and Proposition 15.41 is directly applicable 
and gives the result. 



Proof of Proposition 2.4 Put (^{t) = ht — (f){t). Recall that (p'if) obeys the estimate 
fICTD . Choose G M such that < 6/2, Vt > t*. Hence the function C(t) is 

strictly increasing and 6/2 < < 36/2. Moreover, we choose t^, sufficiently large so 
that 

F{t + s)< V2F{t) for < s < 37r6, and F{t) > V2A2, Vt > t,. (5.8) 

Fix £ G N such that C(2vr(£ + 1)) > t*. Put = C{27i{i + k)), k e N. Note that 
7r6 < Tk+i — Tfc < 37r6. For a given A; G N put 

Fi = min F{t), F2 = max F{t). 



One has 



\<P'{t)\dt < - / Wm'{t)dt 



Tk ^ JC{27T(e+k)) 



2A2AS /-^-C^+fc+i) |coss| 
< — , / , ds 



W^l - ^2 J2nil+k) ^(5) + \ F{sy - S(S)2 sins 



where F = F o (-\ a = a o (-K Put M+ = 27r{i + k) + [0,7r], M_ 
2n{i + k) + [n,2n]. One has 

cossl , f'^^'^ coss , 2 log 2 

' ds < 2 / , ds < ^ 



A/+ F(s) + ■\/f(s)2 - a(s)2 sins ^1 + V^i'-^l sins V^i' - ^1 

For s G M_ one can estimate 



2 



1 ^ F2 + v/:f?^ 



sm s 



+ Jf(s)2 - a(s)2 sin s ^ + ^1 ^ ^1 + V^i - Al sin s 



where we have used that 



< , F^ - Al < 2F^ - Al < AF^^ - 5A^ < A{F^^ - A^), 
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as it follows from f lS.Sp . Thus one arrives at the estimates 



cossi , f^^"^ coss 



^- F{s) + \jF{sY - a(s)2 sin s Jo F^- y^F, - sin s 

and 

Referring to the asymptotic behavior fl2.28p one concludes that there exists a constant 
> such that 

k=l -''^k j=£+l 

Hence the limit limt_5.+oo (f>(t) = 0(oo) G M does exist and (12.291) follows. 

6 Some derivations related to the guiding center 
coordinates 



We keep the notation introduced in Subsection 12.41 but for the sake of simplicity we 
again put m = e = 1. The physical constants can readily be reestablished if necessary. 
In particular, X = q + {l/b)v^, R = —{l/b)v-^, and a direct computation yields 

4 
Using fl2.4p . fl2.5p one derives the equalities 

|XP = i (2/ + \a\ - a), = \(2I+ \a\ + a) (6.1) 





and 



r = (|Xp + + 2|X||i?| sui{if)f'^ , Pr = - \X\\R\ cos(¥?). (6.2) 



r 



On the other hand, one has 

2 



X|2 + + 2X ■ i? = - (^2/ + |a| + 2v//(/+ |a|) cos(?9 - x)) ■ (6-3) 



By comparison of (16.31) with (16. 2p one shows equality (I2.34p . 

Further we sketch derivation of the the asymptotic relations (12.401) . To this end, 
let us compute the derivative x'it)- This can be done by differentiating the equality 

r (cos 9, sin 9) = |X| (cos x, sin x) + |_R| (cos ^, sin {}) 
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and then taking the scalar product with the vector (— sin 6', cos^). One has 

^1 dH a b 

where H is the Hamiltonian fl2.ip expressed in polar coordinates. Using fl2.34p and 
some straightforward manipulations one finally arrives at the differential equation 

^ |X|6r2 ^ ^ 

Equation (16.41) admits an asymptotic analysis with the aid of similar methods as 
those used in Section O In order to spare some space we omit the details. We still 
assume that 0/6 e N. Recalling f l2.32p . f l2.33p one observes that the main contribution 
to the growth of xit) over a period T = 27r/6 equals 

> > "^"^ > ACnT 0.^0 J, 1 + v/r^ sin(6t + 0(00)) 

Proceeding this way one finally concludes that the first equality in f l2.40p holds, with 
x(c>o) being a real constant and 



47r/'(-0 7o V \b J \ b Vyi-cos(t) 

Being given the Fourier series f'{t) = YlT=i (flfcCOs(A;t) + bksin{kt)) one can also ex- 
press 

^ CXI I 00 

D = - ^ (flfc sin(A;^) + hk cos(A;^)) / ^ (ofc cos(A;0 - h sin(A;^)) 
fc=i ' fc=i 

as it follows from the equality, valid for any n G N, 

sin(nt) sin(t) 



1 — cos(t) 



dt = TT. 



Moreover, fl2.32p . fl2.34p and the first equality in f l2.40p imply the second equality in 
(I2.40p where one has to put 'd{oo) = 0(oo) + x(oo) — tt/2. 
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